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The ground-state phase diagram of the half-filled two-leg Hubbard ladder with inter-site Coulomb 
repulsions and exchange coupling is studied by using the strong-coupling perturbation theory and 
the weak-coupling bosonization method. Considered here as possible ground states of the ladder 
model are four types of density-wave states with different angular momentum (s-density-wave state, 
p-density-wave state, d-density-wave state, and /-density-wave state) and four types of quantum 
disordered states, i.e., Mott insulating states (S-Mott, D-Mott, S'-Mott, and D'-Mott states, where 
S and D stand for s- and d-wave symmetry). The s-density-wave state, the d-density-wave state, 
and the D-Mott state are also known as the charge-density-wave state, the staggered-flux state, and 
the rung-singlet state, respectively. Strong-coupling approach naturally leads to the Ising model in 
a transverse field as an effective theory for the quantum phase transitions between the staggered- 
flux state and the D-Mott state and between the charge-density-wave state and the S-Mott state, 
where the Ising ordered states correspond to doubly degenerate ground states in the staggered-flux 
or the charge-density-wave state. From the weak-coupling bosonization approach it is shown that 
there are three cases in the quantum phase transitions between a density-wave state and a Mott 
state: the Ising (Z2) criticality, the SU(2)2 criticality, and a first-order transition. The quantum 
phase transitions between Mott states and between density-wave states are found to be the U(l) 
Gaussian criticality. The ground-state phase diagram is determined by integrating perturbative 
renormalization-group equations. It is shown that the S-Mott state and the staggered-flux state 
exist in the region sandwiched by the charge-density-wave phase and the D-Mott phase. The p- 
density-wave state, the S'-Mott state, and the D'-Mott state also appear in the phase diagram 
when the next-nearest-neighbor repulsion is included. The correspondence between Mott states in 
extended Hubbard ladders and spin liquid states in spin ladders is also discussed. 



PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm, 71.30.+h, 74.20.Mn 



I. INTRODUCTION 

Ladder systems have been studied intensively over the 
years as a simplified model system that shows variety of 
quantum phenomena due to strong electron correlations .EJ 
Since the ladder models can be analyzed with powerful 
nonperturbative methods such as bosonization and con- 
formal field theory as well as with large-scale numerical 
calculations, they provide a useful testing ground of var- 
ious theoretical ideas developed for the two-dimensional 
case. Moreover, the studies of ladder systems have 
been strongly stimulated by experimental developments 
in synthesizing compounds with ladder structur|fi n tiiat 
show superconductivity and spin- liquid behavior .oBta A 
good example is the ladder compound .Sri4Cu2404i that 
shows d-wave superconducting ordertj under pressure 
with Ca doping and charge-densitypatave (CDW) order 
as recently suggested experimentallyfllil Theoretical stud- 
ies on dopedJaxidj^mo^gls^irh as t^pTIuhbard and t- 
J laddersMmelailey BmEla B have estab- 
lished that the dominant correlation is indeed a d-wave- 
like superconducting order, a feature that is reminiscent 
of the d-wave superconductivity in high-T c cuprates. On 
the other hand, undoped half-filled Hubbard and Heisen- 
berg ladders are insula tors tha fc4i a^jo a g ap in both charge 
and spin excitations.lallMl3ByyE2 This spin-liquid 
behavior is caused by singlet formation on each rung, and 
the state is said to be in the rung-singlet phase. It is also 



named D-Mott phaseil because of its close connection to 
the d-wave-like paring state. 

Recent theoretical interest on the ladder models has 
been focused on the search of exotic phases in thepe 
systems. In particular, the staggered-flux (SF) sf. 
which is also known as the pc 
and the . dnd e. n S) ty, . w ? ve f? 1 l 



27 



al antiferromagnetl 
J ' has received a lot 



2'!. SO 



of 



attentionJ±30at!3tiJt!a For more than a decade the SF 
state has been intensively studied in connection with 
the pseudDpaaa. rihijfie Ja— the two-dimensional high-T c 
cupratesBSBfeiaa-BEl The SF state has sponta- 
neous currents flowing around plaquettes, breaking the 
time-reversal symmetry. Even though ladders are one- 
dimensional (ID), the long-range order of the SF corre- 
lation is possible at half-filling, since the symmetry bro- 
ken in this state is discrete. This point was emphasized 
recently in Ref. Q where it is also suggested that the SF 
phase should occur in the phsse diagram of the SO (5) 
symmetric Hubbard modelao Besides the SF phase, 
the ground-state phase diagram of the ladder models can 
include-,the D-Mott phase mentioned above, the CDW 
phase Jl3 and other phases. 

Motivated by these developments, in this paper we at- 
tempt systematic exploration of the ground-state phase 
diagram of a generalized two-leg Hubbard ladder at half- 
filling that has not only repulsive on-site and inter-site 
interactions but also antifcrromagnetic (AF) exchange 
interaction and pair hoppings between the legs. To map 
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out the possible phases in the parameter space of the 
model and to analyze various quantum phase transitions, 
we employ both the strong-coupling perturbation theory 
and the weak-coupling bosonization method. We find 
that the inclusion of the additional interactions leads to 
emergence of various new phases. 

In the strong-coupling approach, we describe the SF 
state as an AF ordered state of pseudo-spins that repre- 
sent currents flowing on the rungs. The effective theory 
near the phase boundary between the SF state and the 
D-Mott state is then found to be the ID Ising model in a 
transverse field. The D-Mott phase is thus interpreted as 
a disordered state of the Ising model. We also present a 
similar mapping to the ID quantum Ising model for the 
quantum phase transition between the CDW phase and 
the S-Mott phased Here the CDW state and the S-Mott 
state correspond to the ordered and quantum disordered 
states of the Ising model, respectively. Furthermore, we 
show that a low-energy effective theory near the phase 
transition between the D-Mott and the S-Mott phases is 
the XXZ spin chain in a staggered field, which exhibits a 
U(l) Gaussian criticality. 

In the weak-coupling limit, we follow the standard 
approach of taking continuum limit and bosonizing the 
Hamiltonian. We obtain a coupled sine-Gordon model 
for four bosonic modes (charge/spin & even/odd modes) 
and analyze it by perturbative renormalization-group 
(RG) method and a semiclassical approximation. The 
scaling equations we derive are equivalent-to those ob- 
tained earlier by Lin, Balents, and Fisher.CJ We depart 
here from the earlier work. We consider four types 
density- wave states with different angular momentum 
s-density wave (= CDW), p-density wave (PDW, which 
is equivalent to the spin-Peierls state), d-density wave 
(= SF), and /-density wave (FDW). These density- wave 
states break Z2 symmetry and can have long-range order 
at zero temperature. We find that in general there should 
appear four types of Mott insulating phases (called S- 
Mott, D-Mott, S'-Mott, and D'-Mott states), each of 
which can be obtained as a quantum disordered state 
from one of the four Z2-symmetry-breaking density-wave 
states. We then study quantum phase transitions among 
these 8 phases and show that a transition between a 
density-wave state and a Mott state is either second or- 
der (in-, the Ising or SU(2)2 universality class) or first 
order.cll Phase transitions between density-wave states 
and between Mott states are U(l) Gaussian criticalities. 
After classifying the phases and the quantum phase tran- 
sitions, we determine the ground-state phase diagram of 
the extended Hubbard model with extra inter-site repul- 
sion and the exchange interaction. We find that the S- 
Mott and the SF phases appear in the parameter space of 
couplings where the D-Mott and the CDW phases com- 
pete. We also show that the next-nearest-neighbor re- 
pulsion stabilizes the S'-Mott state and the PDW state; 
the latter state is connected to the D-Mott state through 
the SU(2) 2 criticality. 

This paper is organized as follows. In Sec. [iljthe model 



we analyze in this paper is introduced. In Sec. Ill we 



study the ground-state phase diagram by the strong- 
coupling perturbation theory, and examine phase tran- 
sitions between the competing ground states: the SF, D- 
Mott, CDW, and S-Mott states. In Sec. [[V]we apply the 
weak-coupling bosonization method to study the ground- 
state phase diagram. We derive effective low-energy the- 
ory for the charge mode and for the spin mode that de- 
scribe the Gaussian, Ising, and SU(2)2 criticalities. The 
connection of our results to the phase diagram of spin 
ladders with spin liquid ground states is also discussed. 
We then determine the phase diagram of the generalized 
Hubbard ladder from perturbative RG equations. Fi- 
nally, the results are summarized in Sec. 



II. MODEL 

We consider a half-filled two-leg Hubbard ladder with 
on-site and inter-site Coulomb repulsions and rung ex- 
change interaction. The Hamiltonian we study in this 
paper is given by 

H = Ht {l + H t± + H mt + Hv l{ + Hy + H pail . (2.1) 

The first two terms describe hopping along and between 
the legs, respectively: 

h h = -*uE( c W c i+i.^ +H - c -)> ( 2 - 2 ) 



Hi 



I. '2.,., • II-''-)- (2-3) 



3>a 



where Cj^^ annihilates an electron of spin cr(=T , J,) on 
rung j and leg l(= 1,2). The Hamiltonian H lnt = 
Hjj + Hy ± + Hj ± consists of three terms representing 
interactions within a rung: the on-site repulsion, 



3,1 



(2.4) 



the nearest-neighbor repulsion on a rung, 

H v± = V± E n. hl n ja , (2.5) 
i 

and the nearest-neighbor exchange interaction on a rung, 



Hj ± — J± 22 SjA ■ s. 



3,2- 



(2.6) 



The density operators are rij.i^ = c ]i a- c j l,a- ano - n i>l 
%,/,T + n j,i,li an d the spin-i operator is given by 



s« = \ E i 



l,<Ti <Tct 1.°'2 C j,l,cr 2 



(2.7) 
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where cr aiM2 are the Pauli matrices. The Hamilto- 
nian (2.1) also has nearest-neighbor repulsive interaction 
within a leg, 



Hv \\ = V \\ ^2 n o\l n o+hh 
o,l 



(2.8) 



and next-nearest-neighbor repulsion, 

Hy = V'/, (n-j.i nj+1,2 + ri. h2 ■ ( 2 - 9 ) 



The last component of the Hamiltonian (2.1) is the pair 
hopping between the legs, 

#pair = tpair ^ (4,1,T 4,14 C 02A C 0,2,1 + H ' C -) ■ ( 2 - 10 ) 



The coupling constants, U, V±, Vn, V, Jj_, and i pa ir, 
are assumed to be either zero or positive. (Most of our 
discussions are actually concerned with the case V\\ = 
V = t 



pair 



half-filled case where ; rijj, equals the number of total 
lattice sites. 



III. STRONG-COUPLING APPROACH 

In this section, we perform strong-coupling analysis 
starting from the independent rungs and discuss transi- 
tions between various insulating phases. 

We begin with eigenstates of H int for decoupled rungs 
at half-filling. Convenient basis states for two electrons 
on a single rung (e.g., jth rung) with Sj X + Sj 2 — arc 



Wo = 


T\ 
1/ 


> = c },i,T c kil°). 




(3.1) 


12), = 


l\ 
t/ 


3 


(3.2) 


13}, = 


n 


' 


(3.3) 


14), = 


n 


' 


(3.4) 



The interaction Hamiltonian H m t is diagonalized as 



11)^-12), 
V2 



V ± - -J± 



II), + 12), / 1 
Hint V2 = r ± + I Jj 
ff int |% - C/|3)„ 



- 12) 



V2 
+ |2)i 



J , (3.5) 

, (3.6) 

(3.7) 
(3.8) 



Comparing the eigenvalues, we find that the lowest- 
energy state of Hi nt for U > V± — 3Jj_/4 is 



|D-Mott> = II 7f 



(3.9) 



This state is a direct product of rung singlets and is noth=. 
ing but the strong-coupling limit of the D-Mott phascO 
or the Mott insulating phase of a half-filled Hubbard lad- 
der. 

When U < V± — 3Jj_/4, on the other hand, the dou- 
bly occupied states |3) and |4) become lowest-energy 
states. In this case, one of the possible ground states is 
the ortsite paired insulating state realized in the S-Mott 
phasejH3 



is-Mott) = n 7f 



TI 



(3.10) 



Another possible ground state is the CDW state: 



0.) In this paper we consider only the and 



|CDW)i = XI 



|cdw) 2 = J; 



2j-l 



TI 



2j-l 



Ti- 



ll 



■2j 



(3.11a) 



(3.11b) 



In the next subsections we study phase transitions be- 
tween these phases. 



A. CDW— S-Mott transition: Ising criticality 

In this subsection we jdiscuss the phase transition, be- 
tween the S-Mott phased and the CDW phaseBEl for 
U < V± — 3Jj_/4. This can be analyzed by mapping the 
system onto an effective spin model. A similar analysis 
for the SO (5) symmetric ladder is reported in Refs. ^4] 
and (H 

We restrict ourselves to the lowest-energy states |3) 
and 1 4) and denote them as 



!+>,• = |3>., !-),- = |4) 5 



(3.12) 



to make the connection to a spin model more evident. 
We regard |±) as the pseudo-spin up/down states. In 
this picture, the antiferromagnetic ordering of the spins 
corresponds to the CDW ordering. We will treat the 
single-particle hopping terms H t ^ and H t± as weak per- 
turbations to derive effective Hamiltonian in the Hilbcrt 
space of |+) and |— ). The lowest-order contributions 
come from the second-order processes: 



H^ = H tj 



1 



Eq — -f/int 
1 



H h , (3.13) 
H t± , (3.14) 



Eq — H int 

where Eq = NU with N being the number of rungs. The 
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nonzero matrix elements of H ( 2a ) and H"( 2h ) are given by 

4t 2 



±,T\H {2a) \± 



(±|ij( 2b )|±) J = (±| J ff( 2fc )|T) J = 



U-2V± 



(3.15) 



2t\ 



U-Vx + 3Jx/4' 

(3.16) 



where |s, s')j = \s) j\s')j + i (s, s' = ±). The above Hamil- 
tonian is written in terms of pseudo-spin operators as 



Jj(2a) = 



#(2&) 



2t 



2V± - U 



(3.17) 



2f 2 



U -V±+3J±/4 



^r? + const., (3.18) 



The coupling constants in the quan- 



tum Ising model are modified to 



K 



2ti 



2V± - U 



2Vn - 2V', 



h = 



2t\ 



V±_ - 3J x /4-U 



(3.23) 
(3.24) 



Thus, Hy.. , Hy>, and i? pa i r do not change the Ising uni- 
versality and only affects the coupling constants. Their 
main effect is to move the phase boundary. The Vy and 
ipair interactions favor the Ising ordered phase or the 
CDW phase, while the V interaction is in favor of the 
S-Mott phase. 



B. D-Mott— S-Mott transition: Gaussian criticality 



where rj and rj are Pauli matrices acting on the pseudo- 
spin states: rj[±)j = ±|±)j and rf\±)j — |=f)j ■ Here we 
find that £f( 2a ) favors antiferromagnetic ordering, while 
jj( 2b ) prevents the order. We thus find that the effec- 
tive Hamiltonian for the doubly occupied states ff£g = 

H (2a) + H (2b) 

is given by the one-dimensional quantum 

Ising model, 



Next we discuss the parameter region U rs V± — 3 Jj_/4. 
In this case the low-energy states of H lnt are formed out 
of (II),- - |2) J )/72, |3)„ and |4),; see Eqs. 
The analysis in the previous subsection indicates that, 
among the states made of |3)j and |4)j-, only the S-Mott 
phase can appear for U « V± — 3J±/4 due to the large 
transverse field h. We thus keep only the two states, 



rroff 



E 

3 



(K- 



(3.19) 



|2)j 



V2 



_ |3) J + |4), 
" V2 



(3.25) 



where the antiferromagnetic exchange coupling K and 
the magnitude of the transverse field h are given by 



K = 



2t\ 



2t 2 



2V± - W 



V± - 3J±/4-U' 



(3.20) 



This model exhibits the Ising criticality at K — h be- 
tween the ordered phase (i.e., the CDW phase) for K > h 
and the disordered phase for K < h. The ground state 
in the disordered phase is essentially the eigenstate of 
t x with eigenvalue +1, which is nothing but the S-Mott 
phase: 



It* = +!),- = 



V2 



| S-Mott). (3.21) 



The condition for the CDW phase to appear is given 
in terms of the Hubbard interactions as 



i-(W*u)* 
V±> T=WJW 



4[1-2(W*||) 2 ] 



Jx, (3.22) 



where < t±/t\\ < l/y/2. When t±/t\\ > 1/V2, the 
CDW phase is not realized within our approximation. 

Here we briefly discuss effects of -Hyj, , Hy/, and i? pa ir, 
treating them as small perturbations. The lowest-order 
contributions come from the first-order perturbation, 
= H V{1 + Hy and H^ = H pail , which can 
be written in terms of the pseudo-spin operators as 
H^) = 2V {] E/r? r/ +1 + 1) - 2V' r/ +1 - 1) and 



for each rung and derive an effective low-energy Hamil- 
tonian for these states to study the competition between 
the S-Mott and D-Mott phases. In this basis, H lnt and 
H t± on the jth rung read 



-Hint = 

Ht ± = 



v±-y x o 

u 





-2t 



-2tj 




(3.26) 
(3.27) 



where = '(1,0) and = '(0,1). Since we are 

interested in the region near the level crossing point U — 
V± — 3J_i_/4, we split the Hamiltonian as 



int 



where the unperturbed Hamiltonian ifpg and the pertur- 
bation term H^ s are given by = U J2j ( n j.i,T n jA,l + 
nj,2,t "j,2,j,+ n 2,i n j-2) and H^ s = (Vj_-U)J2j n j,i n j,2 + 
Hj ± + H t± +Hu- Up to second order in H^ s the effective 
Hamiltonian is obtained as i? (0) + + ff (2) : 



(3.28) 



„(o) _ / U 
H i -\ U 



i \ -2t± 



-2tj_ 




H^ = Ht, 



1 



Eq — Hq 



-Ht, 



(3.29) 
(3.30) 
(3.31) 
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where ffW = £\ ff« = £\ fl}^, and E = NU. 

Now we introduce spin- 1/2 operators SJ, Sj, and S? and 
identify the two states 
s tates of the pseudo-sp 
( |3.30 ) is then written as 



V,=3U/2 



and 



with up and down 



s tates of the pseudo-spin SJ. The first-order term i?W 



H 



(i) _ 



(3.32) 



The energy difference between the |±))j states and the 
rung hopping are represented as the longitudinal and 
transverse magnetic fi elds, r espectively. The nonzero ma- 
trix elements of (3.31) are given by 



±,±\HM\±,±))j = 


2fjj 

'IT' 


(3.33) 


±,±\H^\ T ,Th = 


It 2 
+ U ' 


(3.34) 


±,T\H^\±,Th = 


t 2 
1 

2E7" 


(3.35) 


±,T\HW\T,±h = 


t 2 

+ 1_ 
2U 7 


(3.36) 



where \s,s'))j = \s))j\s'))j + i (s,s' = ±). Thus the 
second-order contribution is written in terms of the 
pseudo-spin operators as 



H 



(2) 



•it 



2U 



(3.37) 



From Eqs. fl3.32| ) and ( |3.37| ) we find that, for U « 
Vj_ — 3 Jj_/4, the low-energy effective Hamiltonian iJgg = 
is given by the anisotropic spin chain under 
the longitudinal and transverse magnetic fields: 



rreff _ 



|./ r .v;'.s;;. . - ./" : (,s;'.s" 

- ( hX $j + ^ 



QZ QZ 



(3.38) 



where J x = 5t 2 /U, J yz = M?JU, h x = 4t±, and 
h z = U — V± + 3Jj_/4. We are interested in the case 
where the Zeeman field in the z direction h z is weak. 
When h z — 0, is equivalent to the XXZ model with 
the exchange anisotropy A = J x / J yz — 5/3r-a,nd a uni- 
form field in the z direction. It is knownE§EJ that the 
XXZ model is in the massless phase governed by the 



V,=U 




FIG. 1: Strong-coupling phase diagram of Ht« +H t± +//mt at 
t± = t||/2 and J x = 0. The CDW-S-Mott transition is in the 
Ising universality class, while the S-Mott-D-Mott transition 
is in the U(l) (Gaussian) universality class. The CDW (S- 
Mott) phase corresponds to the ord ered (disordered) phase in 
the effective quantum Ising model (3.19). The S-Mott and D- 
Mott phases are the ferromagnetically ordered phases of the 
effective spin model ( 3.3S ) . 



c = 1 conformal field theory (CFT) with a compactifi- 
cation radius R (1/2^/n < R < if the uniform 

field is in the range 0.175J" 2 < h x < \j yz . The weak 
perturbation h z is acting on this gapless system. From 
the transformation Sj' z — ► (— iy Sj' z we see that the 
Zeeman field h z acts as a staggered transverse field in 
the antiferromagnetic XXZ model. Since the scaling di- 
mension of (—iyS y ' z is irR 2 , it is_a relevant perturbation 
leading to the opening of a gap .Ell 

Hence we find that, when h z ^ 0, the h z term is always 
relevant and generates a mass gap, while for h z = the 
system reduces to the c = 1 CFT or the Gaussian model. 
Therefore the D-Mott-S-Mott transition is a Gaussian 
U(l) criticality with the central charge c = 1. The critical 
point is at h z = 0, i.e., 



U-V x + -J ± 



0. 



(3.39) 



The character of the gapped phases at h z ^ is deduced 
by looking at the dominant /i 2 -term. Since the gapped 
phases should correspond to states minimizing the rele- 



vant /i z -term, — h z ^ ■ S z , in Eq. ( 3.38 ), we conclude that 



for h z > (h z < 0) the ground state is a ferromagnet- 
ically ordered state with positive (negative) magnetiza- 
tion (S z ), or equivalently, in the D-Mott (S-Mo tt) p hase 
in the original Hubbard ladder model; see Eq. ( fc.25| ). 

The phase diagram obtained from the strong-coupling 
perturbation theory is shown in Fig. |l|, where parameters 
are taken as t±_ = tu/2 and J±_ = 0. The phase transition 
between the D-Mott state and the S-Mott state is de- 
scribed as the Gaussian criticality, while the phase tran- 
sition between the S-Mott state and the CDW state is in 
the universality of the Ising phase transition. The phase 



6 



V,/U 




FIG. 2: Strong-coupling phase diagram of Ht,, +H t± +-ffint at 
tx = t\\/2 on the plane of V±/U and J±/U. The CDW_phase 
occupies the parameter region where the condition 
satisfied. 

diagram for nonzero J± is shown in Fig__§. The CDW 
phase is realized when the condition ( [3.22 ) is satisfied. 
We note that, within the strong-coupling expansion to 
second order, the CDW phase does not exist for tn = t±. 

Finally we discuss effects of the remaining interactions, 
Hv^, Hyi, and £f pa ir- We find that we may ignore Hy t 



and Hy since they yield only a constant energy shift in 
the second-order perturbation theory. By contrast, the 
pair-hopping term changes the phase boundary. Since 



H 



pair | 



I j = and ff pair | 



= t 



pair. 



— ))j, the interac- 



tion part of the Hamiltonian Eq. ( 3.26 ) is modified as 

Hint 



# P air, where 



The main e ffect o f t c 
h z in Eq. (^38|) to 



V± 



Ji 





u + t 



pair 



(3.40) 



pair is to change the coupling constant 
h' = U-V X + 3J ± /4 + t paiT . In 
this case, the critical behavior is still governed by the 
Gaussian theory, and the critical point appears at 



U-Vj 



3 

- J± + t pa 



0. 



(3.41) 



Thus, for tp a i r > 0, the pair hopping term tends to stabi- 
lize the D-Mott phase. As shown in the last subsection, 
it also stabilizes the CDW phase, and the net effect of the 
pair hopping is to suppress the S-Mott phase sandwiched 
by the D-Mott and the CDW phases. 



(|2.10D , The eigenstates of -ff pa i r are given by |l)j, \2)j, 
(13), + 14)^/^2, and (|3), - \4)j)/V2, satisfying 



H 



-ff p air|l), 

|3), ~ |4) 3 - 
;,air V2 
|3)j + \*h 



-ffpair|2)j — 0, 



= -t 



V2 



+t 



|3), ~ K), 
|3)j + K), 



V2 



(3.42) 
(3.43) 

(3.44) 



We thus find that the pair hopping term favors the on- 
site singlet state (|3)j — |4) J )/v / 2- Anticipating compe- 
tition between the on-site singlet state and the rung sin- 
glet state - |2) i )/v / 2 that has an energy gain of 
— 3J^/4 from the exchange term Hj ± , we will consider 
in this subsection the situation where i pa i r — 3Jj_/4 and 
J± is the largest energy scale in the problem. Introduc- 



pair 



ing St 
and H' by 

H 
H' 



where -ff p aj r and H pail are obtained from i? pa ir by re- 
placing ipair with 3Jj_/4 and <5i pa i r , respectively. The 
unperturbed Hamiltonian Hq has eigenstates, 



< P air - 3Ji/4 (Impair | < J±), we define H 



H J± +H, 
Hu 



(0) 



pair ' 

Hy, + H t + H t 



H' 

12 pair i 



(3.45) 
(3.46) 



Ho 


V2 




|1), " |2)j 
V2 ' 


(3.47) 


Hq 


|3)j - |4)i 
V2 




|3),' - |4), 
V2 ' 


(3.48) 


Ho 


+ 12), 
V2 


i 


|1>, + I 2 ), 
V2 ' 


(3.49) 


Hq 


|3>j + |4)i 
V2 




|3), + |4), 
V2 ■ 


(3.50) 



We will focus on the degenerate low-energy states — 
|2)j)/v / 2 and (|3) 3 - — |4)j)/V2 and work with the following 
states that break time reversal symmetry, 



IT), 
II), 



J - |2>^) +^(13),- - |4>^)] , (3.51) 
i-|2)i)-i(|3)i-|4)i)]. (3-52) 



We regard them as states with finite current running on 
the jth rung (Fig. ||), as they are eigenstates of the "rung- 
current operator" defined by 



J; 



(3.53) 



SF state as AF ordering of rung-current and 
SF-D-Mott transition 



In this subsection, we study the SF state in the lad- 
der system using the strong-coupling expansion. Our 
starting point is the pair-hopping Hamiltonian i? pa ir 



with eigenvalues ±2, 

J, |T), =+2 |T)„ Jj\i)j = -2 I i)r (3-54) 

We note that J is not a true current operator for Hq due 
to the pair hopping term. 



IT) 



t 



~9~ ~Y 
-6— -<> 



i 



<>- 



FIG. 3: Schematic illustration of the states | |) and | J.). 
The arrow denotes a state with a finite current running in the 
arrow's direction. 



The SF state has a long-range alternating order of | f ) 
and | J, ) or, equivalentrvv. of currents circulating around 
each plaquette (Fig. To verify the existence of the 

SF phase, we derive a low-energy effective theory, in per- 
turbation expansion in H' , for the low-energy states | f }j 
and | which we regard as up and down states of a 
pseudo-spin. In this picture, the antiferromagnetic or- 
dering of the pseudo-spins corresponds to the staggered 
flux phase. The lowest-order contribution in H' comes 
from the nonvanishing matrix elements in the subspace 
of | T ) j and \ 

( T \H'\ T )j = U \H'\ i ) 3 = \{U + V ± - 5t pail ), (3.55) 

(T \H'\ I )j = {I |5'|T )j = -\{U - Vx - 5t pail ), (3.56) 

from which we obtain the first-order effective Hamilto- 
nian 

H SF =~(U-V X - ftpair) £ ^ + C ° nSt '' ( 3 ' 57 ) 



where cr" are the Pauli matrices (a = x, y, z). The lowest- 
order contributions in ill and t± come from the second- 
order processes, 



H, 



H, 



(2a) 
SF 

(26) 



1 



SF 



Eq — H 
1 

Ht ± — ~- Ht ± 

En — Hn 



(3.58) 
(3.59) 



where Eq = — 3 Jj_iV/4 with N being the number of rungs 

in the system. The nonzero matrix elements of ffgw are 
given by 



< t, I |i4 a) | T, I )j = ( i, T \hw\1, T )j = - 3^, (3.60) 
where u)j = (m? v =1\ I)- We can thus write 

At 2 

H W = ^±J2(a*a* +1 -l). (3.61) 



r (2a). 



SF 



3J_L 



On the other hand, the nonzero matrix elements of H, 



(26) 
SF 



FIG. 4: Staggered flux state described as a Neel ordered 
state of the pseudo-spin states, 1 1 ) and | J, }. 



(TI^ 6) IT)i = 

= (l\Hg b) \i) J 



from which we obtain 



44 
3J ± 



H, 



(26) 
SF 



3JT^ 



const. 



(3.62) 



(3.63) 



From Eqs. ( [Tfj] ), ( |3.6l[ ), and ( |3.63| ), we find that the to- 
tal effective Hamiltonian is the Ising chain in a transverse 
field, 



m 



SF 



E 

i 



ha* 



(3.64) 



where the antiferromagnetic exchange coupling K and 
the magnitude of the transverse field h are given by 



K 



At 2 
3Jx : 



3J ± 



(3.65) 



This model exhibits an Ising criticality at if = the 
Neel ordered phase (if > \h\) corresponds to the SF 
phase, while for K < \h\ the system is disordered. The 
disordered ground state for h > K > is continuously 
connected with the ground state at h — > 00, i.e., the 
eigenstate of a x with eigenvalue +1. This state corre- 
sponds to the D-Mott state in the original Hubbard lad- 
der, since 



|5* = +i>* = ^ (I T>i + !!>,-) 



1 

71 



ID 



|2>, 



|D-Mott). (3.66) 



Hence we conclude that the Ising disordered phase cor- 
responds to the D-Mott phase. 

It is interesting to rewrite the transverse magnetic field 
h as 

" 3Jl 



U-V± + -Ji 



(3.67) 



The SF phase is realized when the inequality 
16i 2 



< U - Vx + 



^pair 



(3.68) 



3J_L 

is satisfied (assuming i|| — <j_ = t), where we have to 
keep in mind the assumption that t pa 



3 7 



8 



IV. WEAK-COUPLING APPROACH 

In this section, we study the phase diagram of the 
generalized Hubbard ladder, treating the two-particle 
interactions as weak perturbations. To diagonalize 
the single-particle hopping Hamiltonian, we define the 
Fourier transform, Cj^ a (k± = 0) = (cj,i j0 - + Cj ) 2, ff )/\/2) 
e 7 >(fc_L=7r) = (cj.im - Cj t 2, a )/V2, and c a {k) = 
^2je~ lkj Cj ta (kj_)/^/N, where k = (k,kj_) and the lat- 
tice spacing a is set equal to 1. The kinetic energy term 
then becomes 



Ho — H t „ 



>(fc)4(fcK(fc), (4.1) 



where e(k) = —2t\\ cosfc — t± cosk±. For t± < 2tn, both 
the bonding (k± = 0) and antibonding (fcj_ — tt) en- 
ergy bands are partially filled, and their Fermi points 
are located at k = ±k Fi k ± with kp$ = ^ + 5 and 
fcf,7r = § — 8, where 5 = sin (tj_/2i|i). At these 
Fermi points the Fermi velocity takes the common value 
v F = 2*|| [1 - (t_ L /2t\ l ) 2 } 1 / 2 . In the following analy- 
sis we restrict ourselves to the isotropic hopping case 



A. Order parameters 

Let us first define order parameters characterizing in- 
sulating phases studied in this section. We consider the 
CDW, SF, p-density-wave (PDW), and /-density-wave 
(FDW) states as possible density-wave ordered states. 
Their order parameters are written as 



° A = ^E/A(fc)4(fcK(fc + Q) 

k.a 

^Et- 1 )^), (4.2) 



where the operator Jpj denotes a current circulating 
around a plaquette: 



Jp.j = ^E (4.1.(7 C j.2,<? + 4.2,(7 C i+l,2,cr 



C ] + 1.2,cr C j+l,\,a + C ]+l,l,a C j,l,<? 



H.c). (4.5) 



The PDW phase is a Peierls dimerized state along the leg 
direction with inter-leg phase difference 7r, characterized 
by the order parameter, 



CPDW = T E ( C }+1,1,<7 C i,l,v ~ C ) + l-2,a C J,2,<7 + Hc 



(4.6) 

The FDW state is a different kind of staggered current 
states. Its order parameter is 



O 



1 



FDW 



(4.7) 



where the operators J±.j represent currents flowing along 
the diagonal directions of plaquettes: 

J+<j = ' E ( C i+l,2,<r c j,l,CT — c j,l,tr c j+l t 2,trj ' (4-8) 
a 

J-* = 1 E ( C )+l,l,a C j,2,a ~ c ),2.a c j+l,l,a) ■ (4-9) 
a 

The long-range order of staggered currents flowing along 
diagonals of the plaquettes has been examined in a spin- 
less ladder systemH3 

We also introduce order parameters of the s-wave and 
d-wave superconductivity, 



° A = ^E/A(fc)c T (fc)c i (-fc) 



(4.10) 



where A = SCs and SCd, and /sc s = 1 and /scd 
cos k — cos k i . 



B. Bosonization 



where Q = (tt, tt) and A = CDW, SF, PDW, FDW. The 
form factor are given by /cdw = lj /sf = cos A; — 

cosfcj_, /pdw = sinfc, and /fdw = sinfc cosfc^. Order 
parameters for the spin density waves are not considered, 
since their correlations decay exponentially in the bulk of 
the phase diagram of our model. It is clear that the CDW 
order parameter, 



Ccdw = \{ n j,i ~ n J' 2 ^ 



(4.3) 



has no nvanishing average in the CDW states ( 3.11a ) and 
( 3.11b] ). The order parameter of the SF state is 



(4.4) 



We bosonize the Hubbard ladder Hamiltonian in this 
subsection. Following the standard bosonization scheme, 
we linearize the energy bands around the Fermi points. 
The linearized kinetic energy is given by 

Ho = VF{pk-k F ^)cl a {k)c p ^{k), (4.11) 

h,p : a 

where the index p = +/— denotes the right /left- moving 
electron. We introduce field operators of the right- and 
left-going electrons defined by 

= 4f E etkXc pAk, 0), (4.12a) 
VL k 

= -J=Xy fcx Cp,a(fc,7r), (4T2b) 
V L , 
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where L is the length of the system: L = Na. The 
linearized kinetic energy now reads 



Ha 



vf / dx 



"X 



-ip-, k F . kj 

dx 



(4.13) 

where k±_ = (n) for £ = + (— ). 

The interactions among low-energy excitations near 
the Fermi points, Hj = H mt + Hy» + Hy 
written as Hj = J dxTLi , where 



-f^pair ; are 



P-<y Ci=± 



9l\\ ^p,a,Ci ^-P,°X2 1>-P,',Ca 



+ 92± ^t,a,Ci ^-P,-<t,C2 ^-P,~cr,iA <Pp,v,C S 
+ 53|| i>t,a,Ci ^L,C2 ^-P, CT ,C4 <P-p,a,C3 

+ 93± ^laXi ^I-<r,C2 ^-P,-°,U 1>-p,v,<3 

(4.14) 

Here e = (1(3 and e = CiC2- The primed summa- 
tion over Q (i — 1, ... ,4) is taken under the condition 
C1C2C3C4 = +1) which comes from the momentum con- 
servation condition in the transverse direction. The cou- 
pling constants and gf± are related to the original 



coupling constants in the Hamiltonian (2.1): 



J dxHo, where 



«° = j£ H H 



dx 



(4.17) 



Here the suffices p and a refer to the charge and spin 
sectors and r = ± refer to the even and odd sectors. 
The operator H ur (x) is a canonically conjugate vari- 
able to 4> vr (x) and satisfies [<j> ur (x),IL u >r'{%')] = i5{x — 
x') 8 V v i 8 rr i. We then introduce chiral bosonic fields 



4>tr(x) 



4>vr{x)^TT I dx'H ur (x') 
J — 00 



(4.18) 



which satisfy the commutation relations 

U+ 



~( x )> < l > v'r'( x ')] = ±*( 7r /4) sgn(z - x') 5 u y 5 r y and 
„(x),(f)~, r ,(x')] = I '4) S„y 8 r y . The right-moving 
and left-moving chiral fields (/) + (x,t) and (j)~(x,r) are 
functions of r — i{x/vf) and r + OxJvf), respectively, 
where r is imaginary time. The kinetic-energy density 
can also be written as 



T~io — 



E E E 

p— ± v—p^a r— ± 



(4.19) 



We also introduce the field 6 vr defined by 9 vr — — 
4>~ r . The 8 field satisfies the commutation relation 
[<j) vr (x),Q v i r i(x')] = —iTrQ(—x + x')S rtr t, where Q(x) is 
the Heaviside step function. 

To express the electron fields in terms of the bosons, 
we define a new set of chiral bosonic fields 



9 if I 

^ = l t V± + j-J±_ + m iie V[| + ^m ii£ y', 



a 4 
= 17 + l e V± 

a 



(4.15) 



: J_L + Z^pair + n^,eVj| + hlTli^V' 

(4.16) 



with the numerical factors defined by l± — ±1, l±,+ = 
T3, l±.~ = ±1. mi j+ = m3 i+ = —1, 7711, - = W3,- = —2, 
777 2 ,+ = +2, ma,- = +1. We have neglected the so- 
called (74 terms describing the forward scattering pro- 
cesses within the same branch (left-/right-mover), since 
including these terms would only cause nonuniversal 
quantitative dif feren ces to the ground state phase dia- 
gram. In Eqs. ( 4.15 ) and ( 4.16 ), we have estimated the 
coupling constants in lowest order in the interaction of 
the Hubbard model. The higher-order contributions can 
play a crucial role of changing topology of a phase dia- 
gram, if different kinds of quantum criticalities acciden- 
tally occur simultaneously when lowest-order coupling 
constants are used, as is the .-case in the ID extended 
Hubbard model at half-fillingH This is not the case in 
the ladder model of our interest, and we will use the low- 



est order form, Eqs. ( |4. 15| ) and ( f4.16| ). 

We apply the Abelian bosonization methocE3 : E2lEil and 
rewrite the kinetic energy in terms of bosonic fields: Hq 



(4.20) 



where p = ±, s = ±, and £ = ±. The chiral bosons 
obey the commutation relations [tp PjS x(x), ip PjS iy(x')] 
= ip-KSgn(x - x')8 s y8c,,c and [f+,s,C><P-,8>,C'] = 



ilt5 s y 5^y. 



The field operators of the right- and left-moving elec- 
trons are then written as 



^p,°X 



exp(ipkF,k±x + ip<Pp,s,(;), (4-21) 



where s = + for u =\ and s = — for a ={. The Klein 

factors rjvXi which satisfy {?7a-,C> Wv'X'} = 25 CT ,o-'^C,C'! are 
introduced in order to retain the correct anticommuta- 
tion relation of the field operato rs be tween different spin 
and the band index. From Eq. 
tor is given by 



21) the density opera- 



Pp,a,d x ) ='^l,<yX%^X'- 



hl^ {x) - (4 - 22) 



The Hamiltonian and the order parameters., contain 
only products of the Klein factors such asliJL2 T = 
Vl,+ Vl,+ Vh-Vl,-, h a = i]a,+ V<T-, and h' ( = VhCVlX' 
which satisfy T = —h^h^ = +h' + h'_. Since T 2 = +1, 
h 2 = (h') 2 = —1, the eigenvalues are T = ±1, h = ±i, 
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and h! = ±i. We will adopt the following convention: 
T = +1. h a = i, h'^ = iQ. 

In the bosonized Hamiltonian the phase field <j) p _ ap- 
pears in the form cos(2</> p _+4(5x) with S = sin - 1 (t j_ /2i || ) . 
Since tj_ (= t») is not small, we can safely assume that 
the 5 is Tjelevant and the electrons are not confined in the 
legs.E3t3o In this case the cos(2^>„_ + ASx) terms be- 
come irrelevant. We thus discard them as well as other 
terms with higher-order scaling dimensions. The inter- 
action term Eq. ( |4.14| ) reduces to 



Hi 



12 <5 i 9 *^) { d ^vr) 

h \_9c+^= cos20 p+ cos26>, 



•p- 



cos 2(f) p+ 

lc+,a- COs2(/) p 



■g. 

■9, 
9— 



V- 



cos 2<f> a 
cos 2(f) (7 

■ g c+ — cos 2(f> p+ cos 29, 
s+ cos20 p _ cos2</> CT+ 

+ 9c= s - cos 2#p- cc,s 2^- 
+ g-^ — cos 2#p_ cos 2# CT _ 
+ .g s+jS _ cos20 CT+ cos20 CT _ 

+ cos 20o-+ cos 2# CT _ 



(4.23) 



where the coupling constants for the bilinear terms of the 
density operators are given by 



9 P + = Y,(9n+9tl-9?\\), (4.24a) 

e=± 

9 P - = Y,<9t;+9il-9l\l (4.24b) 

e=± 

ff«H- = X>2if -5 2 + ±-5i|), (4.24c) 

^- = 53 c (flJ-ffa±-ff"|). ( 4 - 24d ) 



and the coupling constants for the nonlinear terms are 
given by 



(4.25a) 
(4.25b) 

(4.25c) 
(4.25d) 
(4.25e) 
(4.25f) 
(4.25g) 
(4.25h) 
(4.25i) 



We note that the umklapp scattering (the 53 terms) gen- 
erates cosine potentials that lock the 4>p+ field. 



9c+ — = 


~9 3 ±, 




9c+,s+ = 




+ S 8 || 


9c+,s- = 






9c+ — = 


+9~x, 




9t=, s + = 


~9i±, 




9e=,b- = 






97=— = 






9s+,s- = 


+.?l\ + , 




9s+.— = 


+5iT- 





The coupling constants in Eq. ( 4.23| ) are not inde- 
pendent parameters. Imposing the global spin-rotation 
SU(2) symmetry on the interaction terms Eq. (4.14), we 
find that the relations 



n ++ 
#2|| 


- o ++ 

92± 




+ 9tt 


= o, 


(4.26a) 


9 n 


- 9a± 


- 9lf 


+ SiT 


-0, 


(4.26b) 


g n 


- 92± 


-*tf 


+ 9il 


-0, 


(4.26c) 


g n 


~92± 


~ 5i|| + 


+ 9il 


= 0, 


(4.26d) 


9t\\ 


- 93 f 


-9ll 


+ 9H 


= 0, 


(4.26e) 



must hold. In terms of the coupling constants in Eq. 
(4.23), these relations read 



9<?+ + 9a- 
9<T+ - 9<7- 

9t=,s+ ~ $~,s- 

9c+,s+ — 9c+,s- 



2 9s+,s- 

2 9s+ — 
~9t=J= 

~ 9c+,~s = 



0, 
0, 
0, 
0. 



(4.27a) 
(4.27b) 
(4.27c) 
(4.27d) 



We have ignored Eq. (4.26c) which is the constraint on 
the irrelevant cosine term oc cos(2<^ p _ + A8x). Since the 
SU(2 ) symmetry of the original Hubbard Hamiltonian 
(2.1) cannot be broken, the coupling constants in Eq. 
( 4. 2 3| ) must satisfy Eqs. ( |4.27a| )-( p37d| ) in the course of 
renormalization. 

Finally, the order parameters are written in terms of 
the phase fields: 



Ocdw oc cos 4> p + sin 9 P - cos 4> c 
— sin </>„+ cos 6 - sin 



cos v a - 

■p_(_ cus (7 p _ sin >pa+ sin Q u 



Osf OC COS 0p+ COS # p _ COS 



sin < 



O 



sin I 



pdw oc cos <p p+ cos bp 
+ sin (j> p+ sin 9 p 
Ofdw oc cos 4> p + sin 6* p _ sin <j> a 
— sin <p p+ cos # p _ cos 



COS 9 a - 

sm(p a+ sin0 CT _, 
sin^v + sm9cr- 

COS(f> a + cos 9, 

sin 9„- 



COS ( 



Oscd oc e 



COStf„_ COS( 



1 e 



'p- 
9p+ sin ( 



sin i 



'0+ 



Oscs oc e l0 "+ cos 9 
— ie 



p 

p+ sin 



sin ( 



smc 



cos < 



sin < 



COS ( 



(4.28a) 
(4.28b) 
(4.28c) 
(4.28d) 

(4.28e) 

(4.28f) 



C. Critical properties in the charge and spin modes 

In this subsection, we study the ground state phase 
diagram thro ugh q ualitative analysis of the bosonized 
Hamiltonian ( 4.23 ). First we classify the phases that 
can appear at half-filling, and then discuss (a) the Gaus- 
sian criticality in the charge sector and (b) the Ising and 
SU(2)2 criticalities in the spin sector. 
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1. Classification of phases 

In general all the modes become massive in the ex- 
tended Hubbard ladder at hal f-filli ng. This means that 
in the bosonized Hamiltonian (4.23) cosine terms are rel- 
evant at low energies and that the bosonic phase fields 
are locked at some fixed values (integer multiples o£.7r/2) 
where the relevant cosine potentials are minimized.c3 The 
locked phase fields can be treated as classical variables, 
and the average value of an order para meter is found by 
substituting the locked phases into Eq. (4.28). A nonva- 
nishing order parameter signals which phase is realized. 
We can reverse the logic and find the configuration of the 
locked phase fields for each insulating phase by imposing 
its order parameter to have its maximum modulus. This 
is what we do in the following analysis. 

In the SF, CDW, PDW, and FDW phases the ground 
state breaks a Z2 symmetry. Therefore the order param- 
eter of these phases can have a nonvanishing value at 
zero temperature even in one dimension. In each phase 
the bosonic fields <p p +, P -, 4>a+, and 9 a - are pinned at 
a point where the modulus of the c orres ponding order 
parameter is maximized. From Eq. (4.28) we can easily 
find at which values the bosonic fields are locked for the 
four phases. The result is summarized in Table [j| 

Once the configuration of locked phase fields is un- 
derstood for the SF and the CDW phases, we can also 
find that for the D-Mott and the S-Mott phases using 
the following arguments. On the one hand, we know 
from the strong-coupling analysis that these two insulat- 
ing phases are Ising disordered phases of the SF and the 
CDW phases, respectively, where the 9 a - field is locked. 



On the other hand, the Hamiltonian ( 4.23 ) has some co- 
sine potentials that can lock the (f) a - field. Since the <f> a - 
field is a conjugate field to 9 a -, these two fields cannot 
be locked at the same time. In fact, it is knownO that an 
Ising phase transition must be associated with switching 
of phase locking from one bosonic field to its conjugate 
field. We can thus obtain the D-Mott and the S-Mott 
phases from the SF and the CDW phases by exchanging 
the role of the <p a - field and the 9 a - field, arriving at 
the phase locking pattern shown in Table |. A brief com- 
ment on the connection to the superconducting states is 
in order here. If we ignore the p+ mode for the moment, 
the order parameter of the d-wave (s-wave) superconduc- 
tivity takes nonzero amplitude when the locked phases 
((0 p _), ((f> a+ ), and (0O--)) of the D-Mott (S-Mott) phase 



hia4s- 



are substituted _Ll 
previous resultsElE 1 ' 1 ' 12 
ing, the D-Mott state turns into the 



consistent with the 
n \ that, upon dop- 
d-wave supercon- 
ducting state in the t-J or Hubbard ladder. The effect 
of carrier doping is to make the umklapp term irrele- 
vant and to leave the (j) p + field unlocked. The operator 
g20 P + representing the superconducting correlation then 
becomes quasi-long-range ordered. 

It is possible to construct a disorder parameter that 
characterizes the Ising transitions and that has a nonva- 
nishing expectation value in the D-Mott and the S-Mott 



phases. A candidate operator for the disorder parameter 
is 



fij 

Xi 




(4.29) 



In the weak-coupling limit we take the continuum limit 
and express the operator ( 4.2E ) in terms of the bosonic 
fields. We then obtain 



fij = exp^o-O')]- 



(4.30) 



Indeed, the disorder parameter \Xj takes a nonzero value 
in the D-Mott and the S-Mott phases where the <p a - field 
is locked. In the strong-coupling limit studied in Sec. Ill, 
we may impose the condition that n^i + n,- ( 2 = 2 and 
Sf x + S? 2 = on every rung. Under this condition we 
find that exp(i-|X,) = 1 — \Xf and fij reduces to 



f'j 



\f t t 
i=i L 



(4.31) 



which acts on the pseudo-spin states defined in Sees. 
IIIA and IIIC as — and fJ.j\1)i = for 

i < j. This means that we can write \ij — Yli if and 

= HI d i near the CDW-S-Mott and the SF-D-Mott 
transitions, respectively. They are indeed the disorder 
parameter of the quantum Ising modeEj that describes 
the CDW-S-Mott and the SF-D-Mott Ising transitions. 

Since the PDW and the FDW phases break Z2 sym- 
metry, we can naturally expect that these two phases 
should also have their own Ising disordered phases. We 
shall call them S'-Mott and D'-Mott phases for the rea- 
son that will become clear below. The configuration of 
phase locking in the S'-Mott and D'-Mott phases can be 
obtained from that of the PDW and FDW phases by 
exchanging (0a-) and (9<r-)\ see Table §. We see im- 
mediately that the phase-locking pattern of the S'-Mott 
(D'-Mott) state differs from that of the S-Mott (D-Mott) 
only in the locking of the <j)p+ held shifted by 7r/2. This 
implies that the phase transition between S'-Mott (D'- 
Mott) state and the S-Mott (D-Mott) state is a Gaus- 



sian transition in the 



mode, and that the S'-Mott 



(D'-Mott) state should evolve into the s-wave (<i-wave) 
superconducting state upon carrier doping as in the S- 
Mott (D-Mott) state. 

The nature of the S'-Mott state can be de duced 
through its similarity to the S-Mott state (3.10). We 
first note that, as mentioned above, the S'-Mott state 
is related to the S-Mott state by a ir/2 shift of the 4> p+ 
mode, which is equivalent to translation by half unit cell, 
in such a way that the PDW state is related to the CDW 
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TABLE I: Pattern of phase locking. The * symbol indicates that a bosonic field is not locked. IiS are integers. 



Phase 


(4> P +) 


<<?p-> 


(</W) 




(6» CT -> 


CDW 


f lo +7T/l 


f (/o + 1) +vrJ 2 


f Jo + 7r/ 3 


* 


f /() + 7T/4 


SF 




f it) + 7T/2 


f I + 7T/ 3 


* 


f 7o + 7T/4 


PDW 


-(In + 1) + nL 


-f/n + 1) + 7T/9 


ZL/n + „-/, 




-In + 7T 1 4 


FDW 


f(/ + l)+7T/l 




f /O + 7T/ 3 


* 


f /0 + 7T/4 


S-Mott 


f/o+vr/i 


|(/ + l)+7rJ 2 


f /0 + 7T/ 3 


§/o + nh 


* 


D-Mott 


f /0 +7rJl 


J Jo + 7T ^2 


f /0 +7T/ 3 


§/o + vr/ 4 


* 


S'-Mott 


|(J„ + l)+7r/i 


f(/ + l)+7rJ 2 


f /(J + 7T/ 3 


■|/o + 7T 74 


* 


D'-Mott 


f(/0 + l)+7Tll 


f/o+7r7 2 


§/ +7r/ 3 


f /0 + 7T/ 4 


* 



state. This suggests that the center of mass of a sin- 
glet in the S'-Mott state should be located at a center 
of a plaquette. Noting that cos k cos k± is positive (s- 
wave like) at all the Fermi points, fc = (±(-§ + 5),0) 
and (±(f — <5),7r), of the ladder model, we speculate 
that the singlet-pair wave function (or the symmetry of 
a Cooper pair in the s-wave superconducting state real- 
ized upon doping) is of the form cos k cos fcj_c|(fc)c| (— k) 
in momentum space. In real space this corresponds to 
a linear combination of two singlets formed between di- 
agonal sites of a plaquette. From these consideration we 
come to propose the following wave function as a repre- 
sentative of the S'-Mott state: 



|S'-Mott) = J] J(ct 1)T ct +1 



2,1 S,l,lS'+l,2,T 



+ 42,T c I+i,ia- c ki c I+i,i,T)l°)-( 4 - 32 ) 

This state mostly consists of singlets along the diagonal 
direction of plaquettes but also contains resonating sin- 
glets that are formed by two spins on different legs that 
can be separated far away. 

The D'-Mott state consists of singlets that would turn 
into d-wave Cooper pairs upon doping. Since the singlet- 
pair wave function in the D-Mott state is cos k± in mo- 
mentum space, we expect that the singlet pairs in the D'- 
Mott state should be of the form cos k. In real space this 
corresponds to a linear combination of singlets formed 
in the leg direction. This leads to the following wave 
function 



| D'-Mott) = Yl 



E 

1=1,2 



c j,',T l j+i,',1 



n t c t 
-j,U c j+i.U 



|0) 
14.33) 

as a representative of the D'-Mott state. It is easy to see 
by expanding the product that this state is a resonating 
valence bond state in which some singlets can be formed 
out of two spins that are separated arbitrary far away 
along a leg. However, amplitude of the states having 
such a long-distance singlet is exponentially suppressed 
with the distance between the two spins. 

It is interesting to note that the wave function ( 4.32) ) 
can be constructed from the S-Mott wave function 



( |3-10D by replacing c] , CT 
for I = 1 (2) such that 



with c! , , T , 

■t J 



where I = 2 (1) 



)/y/2. This rule can also be used to construct 



C j.l^ C j+l,Ll 

the wave function of the D'-Mott state (4.33) from that 
of the D-Mott state (|3~9|) . 

Since the d>„_ field is locked in the S'-Mott and D'- 



Mott phases, the operator (|4.30| ) also serves as the dis- 
order parameter in the PDW-S'-Mott and the FDW-D'- 
Mott transitions of the I sing universality class. In fact, 
the disorder parameter (4.30) takes a nonzero value in 
any of the Mott phases and vanishes otherwise. 

The various insulating phases and phase transitions 
among them are schematically shown in Fig. ||. In this 
figure phase transitions between a phase in the left col- 
umn and another in the right column, such as transi- 
tions between the Mott phases, are the c = 1 Gaussian 
criticality. It would be interesting to find an order pa- 
rameter that can distinguish different Mott phases. The 
transitions in the vertical direction within a column are, 
if continuous, either the c — 1/2 Ising criticality or the 
c = 3/2 SU(2)2 criticality. The latter may be replaced by 
a first-order transition. We will discuss these transitions 
in more detail in the following subsubsections. 

A brief comment on the related earlier works is in order 
here. The top four phases (SF, CDW, S-Mott, and D- 
Mott) in Fig. || and the Gaussian and Ising transitions be- 
tween these phases have been found in the weak-coupling 
RG analysis of the SO(5X symmetric ladder model by 
Lin, Balents, and Fisher J23 The misidentification of the 
SF phase with the PDW phase made in this wopk has 
been corrected later by Fjaerestad and Marston.E3 We 
have pointed out the existence of four more phases in the 
generalized Hubbard ladder model and determined the 
universality class of the phase transitions between all the 
8 phases. 



2. Gaussian criticality in the charge degrees of freedom 

First we discuss the Gaussian criticality when all the 
modes except the relative charge mode (p— ) become 
massive at some higher energy scale. This situation is 
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— o — ^5 — cS — ^5 — 
CDW 




the effective theory 



S'-Mott 



D'-Mott 



■* ► Gaussian criticality (c=l) 

- Ising criticality (c=l/2) 

-« ►■ SU(2) 2 criticality (c=3/2) 

or first-order transition 



FIG. 5: Schematic illustration of the phase diagram under 
the global SU(2) symmetry. The phase transitions indicated 
by the solid (dashed) arrows are the c = 1 (c = 1/2) criti- 
cality. The phase transitions indicated by the double arrows 
are either the c = 3/2 SU(2)2 criticality or first order; see dis- 
cussion in Sec. IVC3 and Fig. |l^. The diagonal solid arrows 
denote the Gaussian transitions in the 4> p+ mode. 



relevant for the horizontal transitions in Fig. |B|: SF- 
CDW, D-Mott-S-Mott, PDW-FDW, and S'-Mott-D'- 
Mott transitions. We take the D-Mott-S-Mott phase 
transition as an example. Without loss of generality 
we may assume that the phase variables are locked at 
(<j) p +) — (<f>a+) — {4>o-) — mod 7T. Below the en- 
ergy scale at which the three fields are locked, we can 
replace the cosine terms in the Hamiltonian Eq. (4.23) by 
their average: cos20 p+ — > c p+ = (cos2</> p+ ), cos2(j) cr+ — > 
c a+ = (cos2(f> a+ ), and cos20 (T _ — * c a - = (cos20o._), 
where c p+ , c a+ , and c (T _ are nonuniversal positive con- 
stants that depend on bare interactions. We then have 



H p - = V -^[(d x( l> + p _f + {d x ct > --) 2 ] 



+ 



9c 



cos 20 p -, 



2n 2 a 2 

where the coupling constant g— is given by 



9T= = c p+ 9c+ — + c °+ 9c~=. 



s+ 



(4.34) 



(4.35) 



Since the canonical dimension of cos20 p _ is 1, the gzz 
term is a relevant perturbation and hence the system 
always becomes massive except when g-^z = 0. If g-^rz > 0, 
then the phase field is locked as (0 P -) = tt/2 mod w, 
which corresponds to the S-Mott phase. When g-^z < 0, 
the phase field is locked as (9 P -) = mod n, and the 
ground state in this case turns out to be the D-Mott 
state. The Gaussian criticality with the central charge 
c = 1 is realized at grrz =0. In terms of the original 
Hubbard interactions the coupling constant g-^z is given 
by 



.9— 



C [U-V± + -J x +t 



pair 



+ C>(V l{ -V'), 



where C = c p+ + c a+ + c CT _ and C = 
are nonuniversal positive constants. 
(S-Mott) state appears when U — Vx 



2c p+ - 
Thus, 
+ 3J 



C'(Vj| - V')/C 
shows up at 

U-V\ 



(4.36) 
- 2c a+ - c CT _ 
the D-Mott 

7 4 + toair - 



> (< 0), and the Gaussian criticality 



ljx+t paiT -^(V {l -V')=Q, (4.37) 



which is the same as the phase b ound ary obtained from 
the strong-coupling analysis, Eq. ( |3.41 ), for V\\ = V = 0. 

The SF-CDW phase transition can be analyzed in a 
similar way. We consider a situation where the phase 
variable 6* CT _, instead of is locked at (8 a -) = mod 
7T. In this case we can replace the cosine factor in the 
Hamiltonian as cos20 cr _ — > c^— = (c os 26 a -) > 0. The 
effective theory is given by Eq. ( 4.34 ) with the coupling 
constant grrz = c p+ g c+ — + c a+ g— tS+ + c^zg——. The 
SF (CDW) state is realized for g^z < (> 0), where the 
phase Q p - is locked at (7r/2) mod tt. In terms of the 
original Hubbard interactions, the coupling constant g— 
is given by Eq. (4.36) with C = c p+ + c a+ > and C — 



2c p+ +2c CT+ +3c— . We thus conclude that the SF (CDW) 
state appears for U-Vx + \ Jj_+W _C "(^1l -V')/C > 
(< 0), and th e con dition for the Gaussian criticality is 
given by Eq. fl4.37p . 

The other transitions of the c = 1 Gaussian criticality 
can also be analyzed in the same manner. We note that 
in addition to the Gaussian criticality in the p— mode 
discussed above, there is another Gaussian criticality in 
the p+ mode that govern the SF-FDW, CDW-PDW, 
D-Mott-D'-Mott, and S-Mott-S'-Mott transitions. 
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3. Z-2 x 0(3) symmetry in the spin degrees of freedom and 
the Ising and SU(2)2 criticality 



of the SU(2) constraints fl4.40| ), we rewrite the effective 
Hamiltonian in terms of the Majorana fermions: 



Here we focus on the case where the masses of the 
two charge modes (p±) are larger than those of the spin 
modes (o±). Below the mass scale of the charge modes 
we may regard that the (j)p+ and Op- fields are locked 
by cosine potentials. The effective low-energy theory 
is obtained from Eq. ( 4.23 ) by replacing cos20 p+ and 

and 



cos2# p _ by their average values c 



Tin 



■ = (cos 29 p 

VF 
7T 

9a± 

2tt 2 

+ 9 — 
2ir 2 

9s- 



+ 



p+ = 

9s+ 



(cos 2d> p 



2ir 2 a 2 



cos 2<j) a+ 



+ 



27T 2 a 2 

9s+,s- 

2n 2 a 2 
9 S +,J= 
2ir 2 a 2 



cos 2<t> a 



9— 



2ir 2 a 2 
cos2(/> CT+ cos20 CT _ 

cos2(/) cr+ cos2# cr _, 



cos 29„ 



where the coupling constants g s + , 
by 



(4.38) 

and g— are given 



9s+ 
9s- 



c p+ 9c+,s + 
Cp+ 9c+,s- 



"p — 9c— ,s+' 

9c— ,S— ' 



(4.39a) 
(4.39b) 



-9c+ — + Cp-=9z=— ■ (4.39c) 



The coupling constants in Eq. (4.38) are not completely 
free parameters, since the syst em h as the spin- rotational 
SU(2) symmetry. From Eqs. fl4.27| ) and fl4.39| ), the con- 
straints on the coupling constants read 



9s+ - 

9s+,s- 



9s- - g- s 



0. 



9 



1, 
- 2 {9« + 



9a-)- 



(4.40a) 
(4.40b) 

(4.40c) 



To ap preciate the SU(2) symmetry in the effective theory 
(4.38), we fermionize it by introducing spinless fermion 
fields ipp. r (jp = ± and r = ±): 



tp±, r (x) 



exp [±i 2cj>f r (x) 



(4.41) 



where the index r = +(— ) refers to the total (relative) de- 
grees of freedom of spin mode, and {i] r , r) r i} = 26 r . r i. The 

density operators are given by 'ipp ± i>p,± '■ — dx4%±/' !r - 
We then introduce the Majorana fermions £ n (n = 1 ~ 4) 
by 

= ^= (3 + K) > = 7f fe 4 + • ( 4 - 42 ) 

These fields satisfy the anticommutation relations 
#V)} = S(* ~ 0<W<W- With the help 



- i y - e d x £) - im s e + e 



4 



(€+•€-) 



(C + -C-)d-^,(4.43) 



where we have introduced £ p = (£p,£p,£p) and 



m, 



2ira ' 



97= 



2na 



(4.44) 



Thus the effective theory for the spin sector becomes 
0(3) XZ2 symmetric, i.e., the four Majorana fermions are 
grouped into a singlet £ 4 with mass m s and a triplet 
£ with mass m t . We note that the 0(3) XZ2 symme- 
try also appears in the low-energy effective theory of the 
isotropic Heisenberg ladder.EfrE3 It is known that, when 



m s ,m t 7^ 0, the quartic marginal terms lead 
rcnormalization, m s — > fh s and m t — > fh t , wher 



m t = m t 



9<y+ 
2nvp 

3.9cr- 



mt In 



mt In 



A 
\m t \ 
A 

\m t \ 



In 



A 



(4.45) 
(4.46) 



Here A is a high-energy cutoff. The effective theory then 
reduces to 



dxi-) - im t £+ ■ £- 



.V F 



It immediately follows from Eq. ( f4.47| ) that the Ising crit- 
icality with c = 1/2 emerges as m s — ► 0. On the other 
hand, the critical properties for the 0(3) invariant sector 
(fht — > 0) is known to be described by the SU(2)2 Wess- 
Zumino-pNcwikov-Witten model with the central charge 
c = 3/2£M3 

Let us examine the critical behavior in more detail us- 
ing the scaling equations for the cou pling constants ap- 
pearing in the effective Hamiltonian ( 4.43 ): 



dG t 
dl 

dG s 
dl 

dG a+ 

dl 
dG«_ 

dl 



= G t 



= G 
1 



GtG, 
3 



1 



GaGn 



GtG 



2 G ' + 

G a 4-G„ 



1 



2Gi 



(4.48a) 
(4.48b) 
(4.48c) 
(4.48d) 

where dl = da/ a, Gt = —g s +/2'KVF 1 G s = —{g s - — 
gj=)/2irvF, and G a ± = g a ±j2-KVF- The coupling 
G s and Gt are relevant, while G a ± are marginal. 
Within the one-loop RG we find 4 stable fixed 
points, {G* t ,G*,G* a+ ,G* a _) = (±00, ±00,00,00) and 



2GtG s , 
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TABLE II: Signs of the fixed-point coupling constants and 
the masses (m g , fh s , fht) in various phases. 



Phase 




* * * * \ 
,9s-, 9— 9a+,9a-) 


m g 




m t 


CDW 


(+ 


- o,-,+,-) 


+ 


- 


+ 


SF 




-,o,-,+,-) 






+ 


PDW 




+,0,+,+,-) 




+ 




FDW 


(+ 


+,0,+,+,-) 


+ 


+ 




S-Mott 


(+ 


- -,0,+,+) 


+ 


+ 


+ 


D-Mott 




-,-,0,+,+) 




+ 


+ 


S'-Mott 




+,+,0,+,+) 








D'-Mott 


(+ 


+,+,0,+,+) 


+ 







(±oo, T 00 ! °o, ~ °°), which correspond to the 8 phases 
listed in Fig. |^ and Table ||. The Ising criticality is gov- 
erned by the unstable fixed point (G* t ,G* s ,G* <y+ ,G* rJ _) = 
(±oo,0,oo,0), where the Majorana fermion £ 4 is mass- 



less. 



The unstable fixed point (G;,G*,G* + ,G* 



(0, ±oo,0,0) corresponds to the SU(2)2 criticality since 
the triplet £ becomes massless. Finally, we find an- 
other kind of unstable fixed points {G%, G*, G* + , G*_) = 
(0, ±oo, oo, 0), where all the modes are massive. To un- 
derstand the nature of these unstable fixed points, let 
us assume (g s+ ,g— - g s -,g a+ ,g a -) = (0, 2Ai, 2A 2 , 0), 
where Ai j2 are constants (Ai ^ , A 2 > 0). This, to- 
gether with the SU(2) constraint ( 4.40Q , leads to g— = 
—g s - — Ai and g s+ — — g s +,s- — ~A 2 < 0. In this case 



the cosine terms in TL a (4.38) become 

- - ^ o (cos 2</> <T _ - cos26 a -) 

A 2 

- — f T cos20 CT+ (cos20 CT _ +cos26» CT _). (4.49) 
2ir z a z 

Suppose that Ai > a nd (<j> +) = (0 P -) = 0. We then 
find that the potential (|4.49| ) has degenerate minima at, 
e.g., ({</>„+), (4> a _), (0„j) = (0, 0, *) and (§ ,*, §), where 
* means that the phase field is not locked. Since these 
minima correspond to the D-Mott and PDW phases, 
the unstable fixed point describes a first-order tran- 
sition between the D-Mott and PDW phases, respec- 
tively. Hence we conclude that the unstable fixed points 
(Gt,G*,G* a+ ,G* a _) = (0, ±oo,oo,0) correspond to a 
first-order phase transition. The phase transition at 
which the renormalized triplet mass Gjf vanishes can be 
either SU(2) 2 criticality_or first-order transition, depend- 
ing on the sign of G CT +E3. The condition for the SU(2) 2 
criticality is Gt = and G a + < below the energy scale 
where G s becomes of order 1. On the other hand, the 
first-order transition is realized if Gt = and G a + > 0. 

The phase fields are locked at some multiples of tt/2 
depending on signs of the relevant coupling constants at 
a fixed point, (g|—, <7* + , <7|_, <7^—), of the cosine poten- 
tials in Eqs. ( |4.34[ ) and ( 4.38 ). Comparing the config- 
uration of the locked phases and those listed in Table 
H we can find out to which phase the ground state be- 
longs for given combination of the renormalized coupling 



constants, (g±—,g* +: g*_ : g±—). Table || summarizes for 
each phase the signs of these renormalized coupling con- 
stants including g*^, which is positive (negative) when 
(pg-j. (9 a ±) is locked. When writing Table [n|, we have 
used the fact (a) that either one of g*_ and g^— must 
vanish except at the Ising criticality becau se ^ > CT _ and 



Q a - are conjugate fields, and (b) that Eq. ( |4.40a| ) con- 
straints possible combinations of signs of g s+ , g s -, and 



The coupling constants listed in Table |j| also deter- 
mine the signs of masses rn g (= g ^—/2i :a), m s , and fht 
through Eqs. Q4.44Q , fl4.45| ), and ( 4.46Q . The Gaussian 
(c = 1), Ising (c = 1/2), and SU(2) 2 (c = 3/2) critical- 
ities are realized when m a = 0, fh s — 0, and fht = 0, 
respectively. From Table O we can therefore figure out 
which criticality can occur at each phase transition where 
the relevant mass changes sign. The universality class 
of the phase transitions is also summarized in Fig. |^. 
We find from Table that the CDW-S-Mott and SF- 
D-Mott phase transitions are indeed in the Ising univer- 
sality class and the D-Mott-S-Mott phase transition is 
in the Gaussian universality class, in agreement with the 
strong-coupling approach in Sec. III. 

Let us discuss implications of the above general qual- 
itative analysis to the phase diag ram o f the extended 
Hubbard ladder. From Eqs. (4.39) and (4.44) we write 
the bare masses in terms of the coupling constants in the 
model: 



1 

2^ 



2c p+ (U - tpair + V') 

f c^= ( U - V± + j J± + Vir - 4V" 



(4.50) 



U-V x + ^J±+t pail + 2V 



(4.51) 



we assume here that Vji 
is locked at (<£„+) = 



If 



U 



V±_ + 3Jj_/4 > 



To simplify the discussion 
V — tpair = and that 
(mod 7r), i.e., c p+ > 0. 
(< 0), the phase 6» p _ is locked at (tt/2) [see Eq. ( |4.36[ )] 
and c^r = (cos26> p _) > (< 0). Thus, the product 
c^(U — V±+3 Jj_/4) is positive for both positive and neg- 
ative U — V± + 3 Jj_/4, and hence the bare masses m s and 
mt are also positive. We argue, however, that the Ising 
criticality is possible due to the mass renormalization ef- 
fect. The renormalized mass fh s can become negative 
sinc e the coupling constant g CT _ of the correction term in 
Eq. fl4.46|) is given by g„- = 2a(—V± + J±/4). We expect 
that sufficiently large Vi can drive the system toward the 
Ising criticality in the £ 4 mode, even when i pa i r = 0. 

In addition to the Ising criticality at large V±, the 
Gaussian criticality in the # p _ mode should appear at 
V± = U + 3 Jj_/4. Let us find out which phase is realized 
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near the Gaussian critical line. When U — V± + 3 Jj_/4 = 
0, the coupling g a _ equals — 2U — J± and the renormal- 
ized Ising mass becomes 



U 



I- A 



U 
A 



1 



3J ± 

IT 



2Jj 
U 2 



hi 



A 



U + J± 

(4.52) 

where A is a positive constant of order 1. For small 
J±/U this renormalized Ising mass should be positive, 
and we conclude that the D-Mott and the S-Mott phases 
are separated by the Gaussian critical line (Note that 
fh t > 0). As we increase J±/U (or V±/U) along the 
Gaussian critical line, the negative correction (oc g a -) 
in the mass renormalization increases and eventually m s 
can change sign. Across this Ising transition the D-Mott 
and S-Mott phases turn into the SF and CDW phases, re- 
spectively. This implies that a pair of phases surrounding 
the Gaussian critical line changes from (D-Mott, S-Mott) 
to (SF,CDW) at a tetracritical point as J±/U increases. 
This qualitative analysis will be supported in the next 
subsection by a more quantitative renormalization group 
analysis. 

Now we briefly discuss the effect of the pair hopping 
term i pa j r and next-nearest-neighbor repulsion V' . When 
V = 0, the Gaussian transition takes place at U — V± + 
3Jl/4 + tpair = [see Eq. Q4.37| )]. Thus for large i pair , 
we can have a situation where m s < a nd m*_> with 
U-Vx+3J±/4+t pail - [seeEqs. ( fi~50l ) and ([l5l|)l, i.e., 



tp a i r can stabilize the SF state near the Gaussian critical 



line. In the case i pa i r = 0, on the other hand, we expect 
that sufficiently large V can lead to a phase with m s > 
and m t < i.e., the PDW state, if c p+ > c— > 0. 

Finally, we discuss implications of our schematic phase 
diagram (Fig. ||) to the phase diagram of isotropic spin- 
2 ladder systems, which have been studied intensively 
in connection with the so-called Haldane's conjectureO 
about the existence of a finite energy gap in the integer- 
spin Heisenberg chain. By using the abelian bosoniza- 
tion method, it has been shown that four kinds of 
gapped phases can appear in spin laddfit, systems with 
various types of exchange interactionsO'Ed The possible 
gapped phases are (1) the rung singlet state, which is 
known to be realized in the isotropic Heisenberg ladder 
with nearest-neighbor antiferromagnetic exchange cou- 
plings, (2) the Affleck-Kennedy-Lieb-Tasaki(AKLT)-like 
spin liquid state, in which short-range valence bonds cou- 
ple spins on neighboring rungs,Eil (3) the dimerized state 
along chain with it relative phase, and (4) the dimer- 
ized state along chain with zero relative phase. Both the 
rung single state and the AKLT-like state are Haldane- 
type spin liquids with unique ground state and no bro- 
ken local symmetries. In the dimerized states which are 
known to be realized w^?* 1 ! a sufficiently strong four-spin 
interaction is included|a£3 there is spontaneous breaking 
of the translation (Z2) symmetry and the ground state 
is two-fold degenerate. In the limit of large U the ex- 
tended Hubbard ladder we analyze in this paper should 
reduce to a system with only the spin degrees of free- 



dom. This situation corresponds to g-^r < [see Eq. 
( 4.36 )1, i- e -i m g < 0, with \m g \ S> |m s |,|mt|. Under 
this condition, we still have four phases: the SF, D- 
Mott, PDW, and S'-Mott phases. From Table || (see 
also Refs. |54|,[56],[30|) , we can find correspondence between 
the phases in spin ladders and the phases which we have 
obtained in the extended Hubbard ladders: The rung- 
singlet and AKLT-like Haldane states correspond to the 
D-Mott and S'-Mott states, respectively, and the PDW 
(SF) state corresponds to the dimerized state along chain 
with 7r (0) relative phase. We note that physical pictures 
of phases in the extended Hubbard ladder are consis- 
tent with those in spin ladder; for example, the D-Mott 
state is nothing bug the rung singlet state, as seen in 
the strong-coupling approach (see Sec. III). The AKLT- 
like Haldane state, which is known to be realized either 
with plaquette diagonal exchange coupling or with fer- 
romagnetic rung exchange,EJ would be smoothly con- 
nected to the S'-Mott state, in which the ground-state 
wave function consists of single ts fo rmed between diago- 
nal sites of plaquettes [see Eq. ( 4.32| )] and, moreover, has 
the same topological numbers as the AKLT-like Haldane 
state.E3 The PDW state is nothing but the dimerized 
state with interchain phase ir as seen in Fig. |^, which 
is not a Haldane-type spin liquid since the PDW state 
spontaneously breaks translation symmetry and is two- 
fold degenerate. In order to discuss phase transitions in 
spin ladder systems, two kinds of string order parame- 
ters have been introduced which characterize hidden or- 
ders with different topological numbers, i.e., the parity 
of the number aLjdimcrs crossing a line perpendicular to 
the two chainsOEl T hese string order parameters are 
different from fij (Eq. (4.29)), since fij is associated with 
exp(z0 CT _) in the bosonized form while the string order 
parameters introduced in Refs. |6^ and ^32| are associated 
with the cj) a+ field in our notation. Since the phase tran- 
sition associated with the 4>o-+ field is related to fh t — > 0, 
we expect that the string order parameters introduced 
in Refs. |6C] and 62 characterize the SU(2)2 criticahty or 
the first-order phase transition (double arrows in Fig. |^). 
In our schematic phase diagram (|5|) the phase transition 
from the rung singlet state to the AKLT Haldane state 
can take place, (adiich is actually the case in the spin-i 



ladder systemsE2[E3), if the SU(2)2 and the Ising criticali- 
ties appear simultaneously. This implies that the central 
charge for the continuous transition between the rung 
singlet and the AKLT states is given by | + | = 2. This 
transition becomes first order when the marginal interac- 
tion in the triplet Majorana fermion sector is marginally 
relevant. 



D. Renormalization group analysis 

In this subsection, we study the ground-state phase 
diagram of the extended Hubbard ladder model using 
perturbative RG analysis of the 13 coupling constants 
appearing in Eq. ( 4.23| ). These coupling constants are, 
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however, not independent because of th e 4 c onstraints 
coming from the SU(2) symmetry, Eq. ( 4.27 ). Accord- 
ingly, we have 9 independent RG equations that describe 
how the coupling constants scale when we change the 
lattice constant a — ► ae dl . The 9 independent variables 
we choose to work with are: G p + = g p +/2mi F , G p — = 
g p -/2TTv F , G a+ = g a+ /2irv F , G CT _ = g < j-/2nv F , G a = 
(9c+,s- - 9c+—)/ 27TV f, Gp = (5— s _ - g 7= —)/2irv F , 
Ga = g c +—/2nv F , Gb = g c +,s+/2Kv F , and Gc = 
g-^rz s+ /2ttv f . After some algebra we obtain the RG equa- 
tions: 



d G 
d G 

ji Ga+ 

±G 

Ji Ga 



Ji Gb 



Ji Gc 



d G 
dl G0 



+G 2 



A 



2^ 



-G 2 

2 Ga: 



rt2 3 2 1^2 



~2 G *+ ' 
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2 G 
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a — 
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~2^p+ - -jG 



Gb G a + Gc Gp, 
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~2 Gp 
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Gp, 



1 3 
_ 2^p+ G a + ~2~Gb G a 



l G 



G: 
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— Ga Gp, 
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(4.53) 
(4.54) 
(4.55) 
(4.56) 

(4.57) 

(4.58) 

(4.59) 
(4.60) 
(4.61) 



These equations are equivalent to the ones reported in 
Rcf. |2^, in which another set of 9 independent variables 
are used: = (g p+ + g P -)/8, b\ x = -(g„+ + g a -)/2, 



b p 2 



W 4 , 
,p 



'12 



9C, fl2 



(9p+-9p-)/8, fi2 = ~(s<H 
CT _), u p n = -flyi/8, u \2 = 9a/$>, and uf 2 = qb/2, where 
g v = 2ttv f G u . 

Integrating the RG equations ( 4.53| )-( 4.61 ) numeri- 
cally with the initial condition set by the bare coupling 
constants in the extended Hubbard ladder model, we 
find that G p +(l) grows most rapidly and becomes of 
order unity first. At the length scale / = l p + where 
G p +{l p +) — 2, we stop the numerical integration. Be- 
low this energy scale the p+ mode becomes massive. We 
can assume without losing generality that the phase </> p + 
is locked at (<f> p +) = mod tt. The effective theo ry a t 
lower energy scale (I > Z p +) is obtained from Eq. ( 4.23| ) 
through the substitution cos 2(f) p+ — > 1, g c+ — —>■ g^rr, 
9c+,s+ -> 9s+, 9c+,s- g s ~, and g c+ — -> 9^. We 
then derive and solve the RG equations for the coupling- 
constants in the effective theory to understand the low- 
energy properties of the remaining modes. The pattern of 




U/t 



FIG. 6: Weak-coupling phase diagram of Ht,, + H t± + flint 
at t± = tii = t and J± — obtained from the 1-loop RG 
equations. There is a massless mode (C1S0) on the boundary 
between the D-Mott and the S-Mott states while the bound- 
ary between the S-Mott and the CDW state is C0S§. 



phase locking can be found from asymptotic low-energy 
behavior of the g^, g s+ , g s ~, and g— in the numerical 
solution of the RG equations. The phase field $ (= (f> a ± 
or p ^_) is locked at ($) = 7r/2 or 0, if the coupling- 
constant g (g G {g— , g s +,g s -,g—}) behaves as g -> +G 
or — C in the low-energy limit, respectively, where C is 
a positive constant of order unity. Once the configura- 
tion of the locked phase fields is determined, the resulting 
ground state is found from Table [|. The phase diagram 
of the extended Hubbard ladder obtained in this way is 
shown in Figs. |6|-|l0|. We note that this approach repro- 
duces the phase diagram of. the. SO (5) symmetric ladder- 
obtained in earlier studies .E§E3 Since the exotic phases 
like the SF state and the S-Mott state appear only for a 
negative U in this model, we will not further discuss it 
as we concentrate on the case with positive U and V in 
this paper. 

Let us first consider the simple case where U and V± 
are the only electron-electron interactions. The phase 
diagram on the plane of U/t and V±/t is shown in Fig. 
||. In this and other phase diagrams shown below, all 
the modes are gapped everywhere except on the phase 
boundaries. With the standard notation CnSm of rep- 
resenting a state hawig n massless charge modes and m 
massless spin modes,E3 the three_nhases in Fig. || are char- 
acterized as the "C0S0" phase.E3t3 The phase boundary 
between the D-Mott state and the S-Mott state is the 
U(l) Gaussian critical line of th e p— mode (C1S0), which 
is given by V±_ = U; see Eq. ( fj~37j ) with J ± = 0. The 
phase boundary between the S-Mott state and the CDW 
state is the Ising critical line of the spin a— mode, which 



is COS^. This weak-coupling phase diagram is similar to 
Fig. [I] obtained from the strong-coupling approach. 
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FIG. 7: Weak-coupling phase diagram of Ht,, + Ht,+ flint 
at t± = t|| and U/t — 1. This corresponds to Fig. bl Inset 
shows weak-coupling phase diagram of ff t|| + H t± + i/int + 
-Hpair at ij_ = t|| = t, U/t = 1, and t pair /t — 0.5. On the 
boundaries between the D-Mott and the S-Mott states and 
between the SF and the CDW states exists a massless mode 
C1S0. A massless mode C0S| appears on the boundaries 
between the D-Mott and the SF states and between the S- 
Mott and the CDW states. The different choice of U/t does 
not yield qualitative changes to this phase diagram. 



Next, we include the AF exchange coupling J±. The 
phase diagram on the plane of J±/U and V±/U at 
U/t = 1 is shown in Fig. [7| A different choice of U/t does 
not lead to qualitative changes in the J±/U vs V±/U 
phase diagram. An interesting new feature is that the 
SF phase shows up between the D-Mott phase and the 
CDW phase. This is in agreement with the qualitative 
analysis of the previous subsection, where it is found that 
the exchange interaction J± suppresses the S-Mott phase 
and helps the SF phase appear. The Gaussian criticality 
of the p— mode (C1S0) emerges on the almost straight 
phase boundary between the D-Mott phase and the S- 
Mott phase and between the SF phase and the CDW 
phase. This critical line is given by V±/U = 1 + 3J±/AU, 
in accordance with Eq. (4.37). The phase boundary be- 
tween the D-Mott phase and the SF phase and between 
the S-Mott phase and the CDW phase is the Ising criti- 
cality COSi. A tetracritical point of CIS^ appears at the 
point where the two kinds of phase boundaries cross. The 
inset of Fig. ^ shows the phase diagram at i pa i r = 0.5t. 
We see clearly that the pair-hopping favors the SF phase 
over the S-Mott phase. In the strong-coupling perturba- 
tion theory, we have introduced the pair-hopping term 
-ffpair to stabilize the SF state. This is not necessary, 
however, in the weak-coupling approach, where the pair- 
hopping process is effectively generated from the second- 
order process in the rung hopping tj_. In fact, we can 
show that positive pair-hopping terms are generated in 



FIG. 8: Weak-coupling phase diagram of H for U/t = 1, 
Vj| = V± — V , and t pa i r = V' = 0. The tetracritical point 
with C1S| is at (J ± /t,V ± /t) ~ (0.40,0.43). 



the renormalization-group procedure in the SF phase. E3 
Next we turn on the nearest-neighbor Coulomb re- 
pulsion in the leg direction, Vji . The phase diagram 
for Vji = Vj_(= V) is shown in Fig. ||. Even though 
the additional Vji interaction strongly favors the CDW 
state, a small region of the S-Mott phase still remains 
in between the D-Mott phase and the CDW phase. Be- 
sides this quantitative modification the phase diagram is 
not changed qualitatively, and, in particular, the criti- 
cal properties at the phase boundaries are the same as 
in Figs. ^| and Using the density matrix renormaliza- 
tion group method, Vojta et a/.Ej determined the phase 
boundary between the CDW state and a state with ho- 
mogeneous charge density for the model we used for Fig. 
||. At U = 1.5t they observed a transition to the CDW 
state around U/V « 2.9, which is not very different from 
the phase boundary at Jj_ — in Fig. |[ The transition 
is, however, found to be first order for U > At in their 
numerical results, which is different from the continu- 
ous transition we found in the weak-coupling analysis. A 
possible source of this discrepancy might be the neglect 
of irrelevant operators with canonical dimension 4 that 
could become important for strong couplings as in the 
single chain case.EJ 

Finally, we inc lude next-nearest-neighbor Coulomb re- 
pulsion V, Eq. flH|) . Figures | and ^ show the V'-U 
and V-V phase diagrams. In agreement with the discus- 
sion in the previous subsection, the PDW phase appears 
as V is increased. At even larger V the S'-Mott phase 
and the D'-Mott phase appear in Figs. || and |IT^. On the 
phase boundary between the D-Mott state and the PDW 
state appears the SU(2)2 criticality; we have confirmed 
in our numerical calculation that the coupling g a+ in Eq. 
(4.43) is negative, i.e., marginally irrelevant. We have 
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FIG. 9: Weak-coupling phase diagram of H on the plane of 
U/t and V'/t for Vj| = V± = 0, and J± = t pa ,i r = 0. The 
boundary between the D-Mott state and the PDW state is 
C0S§, and the boundary between the PDW state and the 
S'-Mott state is COS^. 



thus established that the two-particle interaction V' can 
drive the system to the SU(2)2 criticality. 

Figure [l^ shows a rich phase diagram containing the 
four Mott phases and the two density- wave phases. We 
note that in Fig. [lO] the six phase boundaries meet at 
V = V = U, which corresponds to C2S2. This happened 
because, with in ou r approximation, all the coupling con- 
stants in Eq. (4.23) except g p + vanish when U = V = V, 



t± = tu, and J± = ip a i r = 0. If t± ^ %ri or if higher-order 
contributions to the t/'s are included,Eil this special situ- 
ation might not occur. In Fig. [l^ the phase boundaries 
between the Mott phases are C1S0 (Gaussian criticality), 
while the CDW-S-Mott and PDW-S'-Mott phase bound- 
aries are COS^ (Ising criticality). The phase boundary 
between the PDW phase and the D-Mott phase is C0S| 
[SU(2) 2 criticality] as in Fig. ^. Finally, the phase tran- 
sition between the CDW phase and the D'-Mott phase 
is found to be first order; we have confirmed that the 
coupling g a+ in Eq. ( 4.43[ ) is positive and marginally rel- 
evant. Even though Fig. [Io| is obtained from the weak- 
coupling RG equations, we think that the phase diagram 
is reliable since we have confirmed that the V/U-V /U 
phase diagram is not changed much when U /t is varied. 



V. CONCLUSIONS 

In this paper we have studied the half-filled general- 
ized Hubbard ladder with the inter-site Coulomb repul- 
sion and the exchange interaction by using the strong- 
coupling perturbation theory and the weak-coupling 
bosonization method. In the strong-coupling approach 
the SF state is described as an AF ordered state of the 



FIG. 10: Weak-coupling phase diagram of H on the plane 
of V'/U and V/U for U/t = 0.5, Vfj = V± = V, and J ± = 
tpair = 0. The phase transition between the CDW and S- 
Mott phases and between the PDW and S'-Mott phases is 
in the Ising universality class (C0S|). The phase transition 
between Mott phases is a Gaussian transition (C1S0). The 
boundary between the D-Mott phase and the PDW phase is 
C0S| [SU(2) 2 criticality]. The transition between the CDW 
phase and the D'-Mott phase shown by the thick solid line is 
a first-order transition. 



Ising model where pseudo-spins represent the currents 
flowing along the rungs. We have shown that the SF 
state can appear next to the CDW state and the D-Mott 
state in the phase diagram and that the quantum phase 
transition between the SF state and the D-Mott state is 
in the Ising universality class. We have also established 
the Ising transition between the S-Mott and the CDW 
phases and the Gaussian transition between the D-Mott 
and the S-Mott phases. In the weak-coupling approach 
we have shown that in general the model can accommo- 
date total of eight insulating phases at half-filling, four 
density- wave phases and four Mott phases (Fig. |J). The 
universality class of the phase transitions among these 
phases is determined. In particular, we have shown that 
the SU(2)2 criticality with the central charge c = 3/2 is 
induced by the next-nearest-neighbor Coulomb repulsion 
V , which drives the system from the D-Mott phase to 
the PDW phase (Figs. | and 0). When V is further in- 
creased, the S'-Mott phase and the D'-Mott phase, which 
correspond to the quantum disordered states of the PDW 
phase and the FDW phase, show up (Fig. ||). 



When this manuscript was almost campleted, we be- 
came aware of the work by Wu et aL,EJ where the 8 in- 
sulating phases in Sec. IV are obtained independently. 
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